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ABSTRACT 

The end state of a coalescing binary of compact objects depends strongly 
on the final total mass M and angular momentum J. Since gravitational 
radiation emission causes a slow evolution of the binary system through 
quasi-circular orbits down to the innermost stable one, in this paper we examine 



the corresponding behavior of the ratio J/M 2 which must be less than l(G/c) or 
about 0.7(G/c) for the formation of a black hole or a neutron star respectively. 
The results show cases for which, at the end of the inspiral phase, these 
C — \ conditions are not satisfied. The inclusion of spin effects leads us to a study of 



precession equations valid also for the calculation of gravitational waveforms. 

PACS number(s): 04.25.Nx, 04.30.Db, 97.80.Fk, 97.60. Jd, 97.60.Lf. 
Subject headings: binaries: close — gravitation 

1. Introduction. 

Binary systems of compact objects like neutron stars or black holes are promising 
sources of gravitational radiation ( [rhorne 1987] ; [Schutz 1996j ); in fact at coalescence they 
emit a great amount of energy in the frequency range of LIGO and VIRGO detectors 
( [Abramovici et al. 1992 ). At least a half of the observable stars come in binary or multiple 



systems ([Batten 1973j ) and a significant fraction of these may evolve to neutron star or black 



hole binaries which would coalesce, due to gravitational wave emission, in a time less than 
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the age of the universe. Moreover, recent estimates ( [rhorne 1987 ; Schutz 1986 ; JN arayan 
Piran, fc Shemi 1991| ; [Phinney 1991| [lutukov fc Yungelson 1993| ; |Yamaoka, Shigeyama, fc 



INomoto 1993| ) make us confident that a few coalescences of compact binaries per year could 
be seen within a distance of about 100 Mpc. 

In order to extract signal from noise, an extremely precise theoretical prediction is 
needed for the gravitational waveform ( |Kidder 1995| ). Therefore much work was recently 
devoted to understand the evolution of the system ( [Kidder 1995| ; pamour 19871 [Lincoln S| 
Will 199Q ; |Kidder, Will, fc Wiseman 1993| ; |Apostolatos et al. 1994| ; |Blanchet, Pamour, &| 
Iyer 1995| , |Blanchet 1996| ; |Will fc Wiseman 1996| ; and references therein) . 



An open question concerns the body coming out of coalescence. It is well known that 
a stationary black hole of mass M and angular momentum J must satisfy the condition 
J/M 2 < l(G/c). Also for a neutron star the ratio cJ/(GM 2 ), hereafter defined as "Kerr 
parameter", is less than about 0.7 ( |Cook, Shapiro, fc Teukolsky 1994| ; |Friedman fc Ipser| 
1992 ; Salgado et al. 1994a| , 1994b), thus it is important to know how this ratio behaves 



during coalescence to see whether a black hole or a neutron star may actually result. 

We also study the precession equations that describe the evolution of the bodies' 
spins. Spin precession plays an important role in the physics of angular momentum, and 
its detailed knowledge is necessary for a consistent calculation of J/M 2 . Moreover spin 
precession is very important for the modulation of the gravitational wave ([Apostolatos et| 
al. 1994| ; |Kidder 1995| ; |Cutler fc Flanagan 1994j ). Our results in this study apply directly to 
this field. 

The coalescence of a binary system consists of two phases. The first phase is a slow 
adiabatic inspiral of the orbits in which the energy and angular momentum loss are assumed 
to be driven only by emission of gravitational waves. The second phase is highly dynamical: 
interactions between the bodies involve their internal structure and full general relativistic 
effects must be taken into account. 

In this paper we consider only the first phase with the purpose of determining the value 
of the Kerr parameter at the end of it, which is the initial value of the Kerr parameter for 
the following phase of evolution in dynamical time. This late-time evolution is not yet well 
understood, because its study is a formidable problem that can be attacked almost only 
by recourse to numerical methods. It is essential for this task to know what the choice of 
initial values can be ( |Cook 1994|) , and this paper is aimed as a help in understanding the 
state of the system in that critical moment of transition between the two phases. 

We adopt the post-Newtonian formalism, which corresponds to a series expansion of 
all quantities in the parameter v/c assumed to be small. There are available in literature 
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formulae for all physical quantities pertaining to the physical system, which are valid 
through second post-Newtonian relative order, i.e., including terms of order {v/c) 4 beyond 
the first nonzero one, and including spin effects ( jKidder 1995| ). This should be a pretty 
good approximation throughout the inspiral phase up to the innermost stable circular orbit 
(ISCO) ( Uunker fc Schafer 1992Q . 



We shall neglect properties of the gravitating bodies connected with their internal 
structure except their spins: no correction terms for tidal interactions will be included 
since they would be equivalent to post-Newtonian corrections of higher order than what 
here wanted ( |Lai fc Wiseman 1996| ). Again this approximation is less justified for compact 



bodies only when their separation becomes comparable to the ISCO radius. 



Our calculations refer to quasi-circular orbits, for it has been proved ( Peters 1964 ; 



Lincoln fc Will 1990|) that all compact binaries have time enough to evolve toward 



quasi-circular orbits before coalescence, except those that are captured at near separations 
in highly eccentric orbits. 

In §|2] we consider the case of nonspinning bodies and examine the behavior of the 
corresponding Kerr parameter. In §[| this analysis is extended to spinning bodies; in 
this case, however, it is essential to take into account the time evolution of the relative 
orientations of the orbital and intrinsic angular momenta. Therefore § [3.1| is focused on this 



topic. Section |37| then gives the results of the analysis of the Kerr parameter when spin 
precession effects are included. Our conclusions are drawn in §|j. We devote an Appendix 
to the study of the precession equations, due to their interest both for the problem at hand 
and more generally for the study of the gravitational waveform, and provide analytical 
solutions for particular cases. 



1.1. Conventions and Units. 

Although we use measurement units in which the constant of light speed in vacuo c 
and the gravitational constant G appear explicitly, we shall always make an effort to work 
with dimensionless physically meaningful variables. 

We refer to a system of coordinates that satisfy the harmonic condition, as usual in 
most related literature, and whose center of mass is fixed by the same "spin supplementary 
condition" (SSC) used in ( Kidder 1995|) , that is the covariant SSC expressed by equation 



(A2a) in that reference. 

We shall denote the "Schwarzschild masses" QDamour 1987|) of the two bodies by m\ 
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and m 2 , their sum by m = mi + m 2 and the system's reduced mass by p = m\m 2 j {rni + mz), 
then set rj = p/m. Following Blanchet, Damour, & Iyer (1995) we define X\ = m\/m, 
X 2 = m 2 /m and notice that, since X x + X 2 = 1 and X\X 2 = i], X\ and X 2 are functions of 
T) only: namely, X x = \±(\- ??) 1/2 , X 2 = \^(\- r]) 1 ' 2 . 

Moreover let x = xi — x 2 be the separation vector between the two bodies in the 
chosen coordinate frame, r = |x|, v = dx./dt, p = (c 2 /G) r/m, 7 = (G/c 2 ) m/r — 1/p, 
Ljv = M x x v )> = Ljv/ILjv'I, L = Ln + post-Newtonian corrections, let E be the 
system's Noetherian energy (Pamour 1987| ), M = m + E/c 2 , and let J be the system's 
Noetherian total angular momentum ( pamour 1987| ), J = |J|. We let S; be the spin angular 
momentum of the zth body, Si = |Sj|, Sj = Si/Si and cr, = (c/G) Sj/m 2 (i = 1,2). We 
write XiX? = &i = \<Ti\, so that Xi is the absolute value of the zth body's spin in units 
of (G/c) m 2 or Kerr parameter, and ranges from to 1 for a (Kerr) black hole and from 
to about 0.7 for a neutron star ( [Friedman fc Ipser 1992 ; Salgado et al. 1994a| , 1994b). 



For purposes of later convenience in writing long equations we also define / = Ljv • si, 
g = Ln ■ s 2 , h = §1 • s 2 , V = Ln ■ §1 x s 2 . 



2. Nonspinning bodies. 



As stated in the introduction, we assume here that the evolution of a close binary 
of compact objects can be approximately described by an adiabatic inspiral of pointlike 
masses through quasi- circular orbits, when the bodies' separation is large compared to their 
sizes ( [Lincoln fc Will 1990| ). The total mass of the system is given for both spinning and 
nonspinning bodies and to the order relevant for our purposes by the following expression 
( Wagoner fc Will 1976t [Junker fc Schafer 1992Q : 



M = m 



l-|7+|(7-r;)T 2 + 0(7 5/2 ) 



In the case of nonspinning bodies (S% = S 2 = 0) the general formula for J ( Kidder 



1995|) reduces to 



J = p (Gmr) 1 ^ 2 Ltv 



l + 2 7 + i(5 - 97/) 7 2 + (7 5 / 2 ) 



Since we have assumed an adiabatic evolution of the system along the sequence 
of quasi-circular orbits, up to the innermost stable one at r ISCO ~ 6Gm/c 2 , the actual 
expressions of the mass and angular momentum loss rate due to gravitational radiation 
emission are not essential to evaluate the behavior of M and J. We can then study directly 
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the behavior of the Kerr parameter. We obtain: 

(2) 



C J 1/2 

GlP =r]p 



1 + (2 + 77)7 + ~ (10 - 17r] + Arj 2 ) 7 2 + 0( 7 5/2 ) 



Since 77 ranges between and |, the latter corresponding to a system of two bodies 
with equal masses, any deviation from equipartition of masses in the binary system would 
decrease the Kerr parameter. 

Figure 1 is a contour plot of cJ / (GM 2 ) in the (p, 7/)-plane. 

The approximation Si = S 2 = is good for the description of systems consisting of 
two slowly spinning neutron stars, since the common values of their Kerr parameters are 
less or much less than 1CT 2 (|de Felice fc Yunqiang 1982|) . We notice that any process of 
slow energy and angular momentum dissipation that keeps the system in quasi-circular 
orbits decreases the Kerr parameter below unity just at the limit of validity of our approach 
(p ~ 6). At this separation the conditions for the formation of a stationary rotating neutron 
star are hardly met; first because J/M 2 < 0.7(G/c) holds only if rj is less than about 0.2, 
which corresponds to a more massive body about 2.6 times as large as the other, and this 
is not the case for a system consisting initially of neutron stars, second because the total 
mass M must be less than the limiting mass of neutron stars. Thus we conclude from an 
inspection of Figure 1 that, since the most probable case for a neutron star binary is that 
they have nearly equal masses, the outcome of a final coalescence could be either a fast 
rotating black hole or a highly rotationally excited neutron star. 



2.1. Estimate of accuracy. 

It is not at all straightforward to estimate the accuracy of equation (Q) for all choices 
of 77 and all values of p. In order to perform a rigorous estimate we should either know 
something about the convergence of the post-Newtonian series (eq. [Q]) or be able to solve 
the problem exactly. Instead our information on the post-Newtonian series is rather poor. 
We neither have an upper bound for the series remainder nor know anything about the 
first neglected terms. In fact we exploited all the information in our hands — that is the 
post-Newtonian terms known up to date — for our calculations. Moreover there are hints of 
slow convergence of post-Newtonian series, reported by some recent studies ( Poisson 1995 



Simone et al. 1997] ) . If we assume that the largest error source in truncating the series is 
given by the first neglected term, we can take the term of order 7 s / 2 in brackets in equation 
(0) as the relative error, or better, as its order of magnitude. Moreover, we shall assume 
that the coefficient of the term 7 5 / 2 is 1, a simplification possible because our units are such 



- 6- 



that all coefficients in the post-Newtonian series are of order 1. Some values for the relative 
error thus estimated are given below: 

error ps 3 • 1CT 8 , at p — 10 3 , 

error ps 10~ 5 , at p = 10 2 , 

error ps 0.11, at p = 6. 

Our confidence in the above estimate may be strengthened by a comparison with the 
case of a test mass in circular orbit around a Schwarzschild black hole. This problem of a test 
mass and a Schwarzschild black hole is well known and exactly solved. Actually it provides 
almost a benchmark for all post-Newtonian calculations, since agreement in the limit of 
rj going to zero gives a necessary condition for their validity. We might also conjecture 
( [Kidder, Will, fc Wiseman 1993|) that this comparison gives a pessimistic estimate of the 



error, since the post-Newtonian series is poorly convergent in the Schwarzschild case. The 
formula analogous to equation (|2|) for a test mass in circular orbit around a black hole reads 



c J 1/2 ^(1 + 7)7(1-27) 



G M 2 



VP 



{l+^(l-7) 2 /[(l + 7)(l-2 7 )]} 



2 ' 



Since both equation (0) and the above equation are when r] is 0, we actually shall divide 
them by r], before any comparison. Thus the second post-Newtonian development of the 
above formula coincides with equation @ at rj = 0. The difference between the two 
formulae — which is a rigorous estimate of the post-Newtonian terms neglected in equation 
(§) in the limit r\ = — is very small, being less that 10~ 3 for p ^ 28.3, and showing a 
sudden increase for decreasing p only at p ^ 16 (due to Schwarzschild behavior) but keeping 
anyway below 6.4 x 10~ 2 for p > 6. 



3. Spinning bodies. 

Most celestial bodies have spin; therefore it is interesting to study how intrinsic rotation 
would affect the evolution of the Kerr parameter for coalescing binaries. 

As in the previous case, we shall confine ourselves to post-Newtonian corrections of 
order 7 2 at most. Since spin-orbit and spin-spin interaction terms will affect the total mass 
only at order equal to or higher than 7 5//2 , inclusion of spin will not change the equation 
(HI) for the total mass. On the other hand, the post-Newtonian corrections to the angular 
momentum will include spin effects already at relative order 7 s / 2 . With the assumption 
of a slow precession of the orbital angular momentum and spins around the total angular 
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momentum, which is valid as long as 7 1 ( Kidder 1995 ), we may take the equation for J 
averaged over one orbit ( |Kidder 1995| ) and rewrite it with our notations as 



Ami J V 4 m 2 



3= fi{Gmr) 1/2 L N \1 + 2 7 - 

+ U(5 - Qiy) - A[(«n . <7 2 ) - 3(L N ■ aJiU ■ <r 2 )]j 7 2 + 0( 7 5/2 )) + 



7 3/2 + 



-Si 



l_I( 3r/ + X 2 )7 + 0(7 2 ) 



+ S 2 



l_I( 3r? + Xl )7 + 0(7 2 ) 



(3) 



We now calculate the square power of J/M 2 from equations @ and ([[]) and obtain 



c J \ 2 _ 



A + A xl x ' 2 + A 27 + A 37 3/2 + A Al 2 + 0( 7 5/2 ) 



(4) 



where: 



r A = v 2 , 

A x = 2r)[L N ■ (en + (T 2 )}, 

2rf{2 + r]) + aj + a 2 + 2(a 1 ■ a 2 ), 
2t 7 [(L jv • o- 1 )(2X 1 + 77) + (Ljv • <r 2 )(2X 2 + V )}, 



A 2 
A 3 
A 4 



(5) 



§77 2 (6 - 377 + 2r7 2 ) - |(1 + r?)^ ■ <r 2 ) - |(1 - 77) (L* • tn)(L N ■ <t 2 ) + 
-^\X 2 u\ + (L N ■ <n) 2 (8 - X 2 )\ - ^[X x o\ + (L N ■ <r 2 ) 2 (8 - Xt)]. 



We shall deal with equation (P rather than try to deduce a post-Newtonian expression 
for J/M 2 as a truncated series because deciding which term in equation (|j) is the leading 
one depends on the choice of the parameters rj, xi and \2 and on the relative orientations 
of Lat, ct\ and cr 2 . As a matter of fact, if we focus on the first few terms of equation ([|), we 
can write 



c J \ 



) = M 1/2 (l + 27)L^v + 0-1 + <r 2 ] 2 + 2r/ 3 + 0( 7 



and notice that the Kerr parameter for spinning bodies depends mainly on the variables 
of the vector r/p 1 / 2 (l + 27)Ljv + <j\ + cr 2 (that is, the first post- Newtonian form of J in 
units of Gm 2 /c) which determine the relative weight of the leading terms in equation (£|). 
For example, if Xi — Xi — 1 an d V = 7 (which corresponds to a more massive body about 
5 times as large as the other), v4 2 7 is of the same order of magnitude as A already at 
values of p ~ 50 and greater at shorter separations. Let us remark that this does not imply 
that the second post-Newtonian formula (eq. 0]) is no longer valid at p ~ 50, because the 
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parameter of development of the post-Newtonian series is 7 which does not depend on 77, 
0\ and cr 2 . We just point out that calculating J/M 2 needs care. 



As for the accuracy of equation (f|), considerations quite similar to those stated in §^T 
about the post-Newtonian series (eq. ||) apply. A thorough discussion should take into 
account the complicated form of the coefficients of the post-Newtonian series (eq. see 
eq. and the unknown similar expressions for the neglected terms). Anyway, since such 
coefficients cannot be much larger than unity, the conclusions we reach are substantially 
equivalent. 

As expected and already noticed, equation depends on the relative orientations 
of the orbital and intrinsic angular momenta of the coalescing bodies, namely, on the 
quantities • 07, L^r • 0*2, and <T\ ■ <x 2 , which are functions of the coordinate time t through 
r. Such functions were considered already by Cutler & Flanagan (1994) because of their 
importance to the secular growth of the gravitational wave phase. We performed a detailed 
analysis of these functions, finding results that we shall report in Appendix, since they have 
more general application than just in our problem. A related and interesting topic is also 
the evolution of the absolute orientations of L^v, 07, er 2 , since the observed gravitational 
waveform depends on the angle between the orbital plane — almost orthogonal to Ln with 
post-Newtonian corrections along cr\ and cr 2 — and the wave propagation direction. This 
question is already fully examined in literature ([Kidder 1995| ; |Apostolatos et al. 1994Q , 
therefore we shall not consider it. 



3.1. Effects of spin precession. 

Here we shall deduce the equations to solve for determining the coefficients of equation 
(|) under the following two assumptions: first, that the spins precess while keeping constant 
absolute values and that the gravitational radiation emission only affects the orbital angular 
momentum L ( [Apostolatos et al. 1994| ; [Kidder 1995| ); second, that the precession frequency 



is much less than the orbital frequency. As analyzed by Kidder (1995) the ratio of the 
precession frequency to the orbital frequency goes as 7. 

We are interested in the dependence on r of the scalar products between any pair of the 
three vectors L N , 07, cr 2 . We put together the equations for the orbital angular momentum, 
separation evolution, and spin precession ( [Kidder 19951 ) averaged over one orbit to lowest 
order and have 

t in M/2f dr 64 3 

L N = niGmr) 1 L N , — = — -0/77 , 

at 5 
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f dL 



N 



G 



dt 



dS\ 
~df 

dS 2 
~df 



2c 2 r 3 



mi J \ m 2 



-3[(L^-S 2 )S 1 + (L iV -S 1 )S 2 ] xL,} 
G 



x L N + 

32 c 3 7] 4 



2c 2 r 3 

G 

2c 2 r 3 



4 + 3— ) L^v + S 2 - 3(Ljv • S 2 )L N 
mi 



4 + 3— ) L w + Si - 3(L W • S^Ljv 
m 2 1 



5 G m 
x Si 



7 L 



A?- 



(6) 



x S,. 



Combining the above equations and rewriting the result in terms of dimensionless 
quantities we obtain: 



~dp • ^ 



15 



d_ 

dp 

d_ 

dp 



L N ■ <7 2 ) 



>i • <r 2 ) 



128?7 

15 
I2877 

15 



1 (Ljv ■ cr 1) 

X~2~ 



1/2 



1 



VP 



(L N ■ CT 2 ) 



VP 



1/2 



[L N ■ <ji x cr 2 ), 



fL^ • CTi X <T 2 ) 



(7) 



Ym 2 


_m x \ 


Ami 


m 2 ) 



12877 



— (Lat ■ <T 2 ) (Lat - CTi X <7 2 



Expanding the quantity [La? x (o"i x cr 2 )] 2 first according to the general rule for 
ax (b x c), then according to that for (a x b) ■ (c x d), we get an equation to close system 

(Lat ■ cti x cr 2 ) 2 = a\al- (cr i- ct 2 ) 2 - al(t N ■ (T 2 f - al(t N ■ (Ti) 2 + 
+2(<Ti • <t 2 )(Lat • cri)(L N ■ a 2 ). 

The system of equations (0) contains nonlinear differential equations whose solutions 
have an oscillatory behavior. In order to avoid that a numerical integration stops when 
Lat ■ <ti x cr 2 first reaches 0, we also differentiate equation (^) and substitute it with the 
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following: 



^-(Ljv • a x x <j 2 ) = {(Xi - X 2 )[(tn ■ a 2 ) - (t N • <ti)(Ljv • o"2)]p 1/2 + 

- (Ljv • <7i) (<Ti ' <X 2 )^-^- + (Ljv • CT 2 )(cTi • <T 2 ) 1 



(9) 



+ (Lat ■ <T 1 )(L A r ■ cr 2 )(L N ■ crx — L N ■ cr 2 ) — — (Ljv • o*i)+ 

A 2 

+-(L W • <r 2 ) + 1. 

We have solved the system of equations (0), (||) numerically for a few choices of the 
parameters xi> X2, X\jX% and of the initial values. Typical behaviors of the solutions to 
the above equations will be shown in the figures that illustrate the Appendix. We have 
judged the analysis of the system of equations (0), (|) explained in Appendix important for 
understanding the role of spin precession during coalescence, since in general the behavior 
of the relative orientations among Ljv, &i, and <r 2 is not simply predictable: for example 
even the average values of the oscillating functions do not keep fixed but evolve, as can be 
seen in Figure 2 — an example in which this effect is dramatic. 



3.2. Results. 

Once the evolution of Ljv • 0*1, Lv • cr 2 ! o~i- cr 2 is known, we can deduce the behavior of 
the Kerr parameter from equations (|j) and (|5|). Plots of this quantity for \i — Xi = 1 an d 
different mass ratios are shown in Figure 3 (solid lines) and compared to the corresponding 
spinless case (dashed lines). 

If the intrinsic spins of the bodies are not very small and if they are directed with a 
positive component along the orbital angular momentum, then the value of J/M 2 remains 
larger than l(G/c) over the entire inspiral phase (Fig. 3a). When one body is much 
more massive than the other and its spin is initially positively oriented with the orbital 
angular momentum, then it contributes to keep the Kerr parameter larger than 1 even if 
the companion's spin is negatively oriented, as shown in Fig. 3b. When both spins are 
negatively oriented with respect to the orbital angular momentum, then the Kerr parameter 
is less than in the corresponding spinless case, at least initially. In fact, at the end of the 
inspiral, the orbital angular momentum has been almost completely radiated away and only 
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the spin is left over, giving a positive contribution to J/M 2 irrespective of its direction. 
For intermediate separations of the two bodies, it may happen that the spin corrections to 
J/M 2 at different post-Newtonian orders add up to 0, as in the case of Figure 3c (notice 
the crossing of the two curves). It may also happen that total spin and orbital angular 
momentum have vanishing sum, as shown in Figure 3d (notice the solid curve approaching 
zero); if this phenomenon occurs before final coalescence, the evolution of the angular 
momenta is known as transitional precession ( [Apostolatos et al. 1994Q . In this last case 



the value of p at which J/M 2 approaches zero varies with the system's mass ratio and 
initial directions of the angular momenta. Thus there are particular choices of these values 
such that the Kerr parameter decreases almost to zero at the end of the inspiral phase; an 
example is shown in Figure 4 (the apparent flattening of the curve at p ~ 6 is uncertain since 
the approximation there reaches its limit of validity). The need for special and uncommon 
choices (see discussion by Apostolatos et al. (1994), §IV.D.l) of the parameters in order to 
obtain a final vanishing total angular momentum makes us think that the formation of a 
Schwarzschild black hole, as result of coalescence of a compact binary, is a rare event. 

The dynamics of inspiralling bodies and of their angular momenta, as described by the 
above formulae, is exactly the same even if the bodies are massive black holes. We can 
infer from our analysis that the end of the inspiral phase of coalescence of two black holes 
in the center of an active galactic nucleus is most likely characterized by conditions which 
are compatible with the formation of a fast rotating black hole. 

Since J/M 2 has a complicated dependence on Ljy Si, Ljv ■ §2, §1 • §2, and the parameters 
f}-, Xij Xii we computed the probability density function of cJ/{GM 2 ) at p = 6 numerically 
(by the Monte Carlo method) for some plausible choices of parameters corresponding to 
neutron star binary systems. The probabilities of mi and m 2 were assumed uniformly 
distributed over an interval ranging from 1.2 to 1.6 (the scale does not matter) and the 
orientations of Ljy, §1, §2 were taken at random, with uniform probability, among all 
directions of space. The parameters \i an d Xi were also taken at random but with a 
probability distribution uniform in logx and ranging from 7.96 x 10~ 5 (corresponding to the 
period of PSR J1951+1123, that is the largest known pulsar period; |Camilo fc Nice 1995|) 



to 0.7 (corresponding to about the highest value compatible with rotational stability). The 
result is the extremely peaked function plotted in Figure 5a. It has a mean equal to 0.86 
with standard deviation 2.94 x 10~ 2 , coefficient of skewness +7.67 and coefficient of excess 
(kurtosis) equal to +6.93. A normal distribution with same mean and deviation is plotted 
(without tails) for comparison (dashed line). 

In Figure 5b the same distribution curve as in Figure 5a is compared to the probability 
distribution of J/M 2 in the spinless case (again numerically calculated by the Monte Carlo 
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method and with the same assumed mass distribution function as before). This latter 
curve, dashed in figure, is even more peaked than the former, and it reaches its maximum 
at the value that equation (Q) takes on for p = 6 and r] = 4, that is J/M 2 m 0.8675(G/c). 

4. Conclusions. 

Neglecting finite size effects we examine the evolution of the Kerr parameter all the way 
up to the state in which the separation between the bodies is 6Gm/c 2 corresponding about 
to the onset of the dynamical instability of circular orbits. We found cases (for example, 
those shown in Figs. 5a and 5b) in which it remains larger than 1 (another instance was also 
found in the numerical study by Wilson, Mathews, & Marronetti 1996). This implies that 
the formation of a stationary black hole is possible in those cases only if J/M 2 is decreased 
efficiently in the final dynamical phases of coalescence. As a matter of fact we cannot infer 
anything on the nature of the coalesced body in such cases. An argument like the one in 
the discussion by Cook (1994), concluding for a final Kerr parameter less than 1 without 
going through the full solution of the evolution problem, simply is not valid: first, because 
it would reach the same conclusion irrespectively of the value of J/M 2 at the beginning of 
the dynamical phase (eq. [16] in that paper always gives J/M 2 < 1); second, because it is 
based on the "a priori" not justified assumption of having a Kerr black hole at the end of 
the coalescence, and this already by itself implies that J/M 2 < 1. 

Similarly the formation of a stationary neutron star seems highly unlikely in the 
absence of a powerfully dissipative mechanism, since we found that in all cases except a 
few very special ones J/M 2 is still greater than 0.7(G/c) at the end of the inspiral phase. 
Anyway, when the final coalesced object is a black hole, our results give us confidence that 
it is rapidly rotating ( |(Jook 19"9"4l ). 

In the course of studying our main object, we had to deal with the problem of 
precessing spins. We found analytical solutions to the precession equations which we 
reported since they have a wider application than just to our purpose: that is, to the 
calculation of the gravitational wave modulation due to spin precession, as fully explained 
by Cutler & Flanagan (1994). 

This work was partially supported by the Italian Space Agency (ASI), the Ministero 
della Universita e della Ricerca Scientifica e Tecnologica (MURST) of Italy and by the 
GNFM of the Italian Consiglio Nazionale delle Ricerche (CNR). We would also thank 
the scientific editor, B. Haisch, and the anonymous referee for suggestions and positive 
criticisms that led to improvements of our paper. 
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A. Analytical solutions of the precession equations. 

A better insight into the behavior of the solutions to the system of differential equations 
(0)' @ i s ensured by an analytical study of them. Because of the complexity of the problem, 
we have been able to obtain only approximate solutions for some cases. We have judged 
this analytical study useful because it applies also to the calculation of the gravitational 
waveform flCutler fc Flanagan 1994]) . 



A.l. Preliminary considerations. 



By elimination of Ljy ■ (Tj x <r 2 from equation (0) we find two relations that do not 
involve L N ■ cr x x CT 2, namely, 



( d d 

p- l ' 2 — {a 1 ■ <x 2 ) + 7]—[L N ■ (o-i + tr 2 )] = 



dp 



dp 



p- l / 2 ^[(t N ■ ai )(t N ■ a 2 )} = ^-[X 2 (L N ■ o-O + X 1 (L N ■ a 2 )). 



(Al) 



Equations (|A1|) are of the form dip/ 'dp + p~ l l 2 dip/ dp = 0, hence an integration by parts 
yields ip + p _1//2 y? + \ f P P~ 3/ ' 2 V 9 (p) dp = 0, and since (p is in all cases at most of the order 
of unity, then for p ^> 1, this last equation leads to "approximate" first integrals in the 
following form: 

ip = constant + 0(p -1 ^ 2 ). 



In order to deal with quantities which have manifestly the same order of magnitude, 
we find it convenient to rewrite the system of equations (0)-(|J) in terms of the variables /, 
g, h, and V: 



dp 

dg_ 
dp 



V 1 



15 (Xl n -l/2\ T/ 

i28{x;- XlX2fp n 



15 (Xx 



128 " XlX29P 



1/2 ) V, 



dh 


15 


(X 2 


--) 


dp 


128 




x 2 J 



Xi 

x 2 



p 1 +Xi—f-X2^r9 



v. 



l-f 2 -g 2 -h 2 + 2fgh. 



(A2) 
(A3) 
(A4) 
(A5) 
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Let us denote the initial values of /, g, h, and p by / , go, h , and p respectively, and 
the initial value of • Si x s 2 , with its sign, by Vq. Moreover we shall denote the initial 
value of the right-hand side of equation by u^Vq. 



A. 2. One spinning body. 

In the very special case where S\ = or S2 = "exactly" , we have 

Ljv • §1 = constant (if S2 = 0) , 

trivially from equation (||). 

A. 3. Equal masses. 

Let Xi = X 2 and thus rj = 1/4. Then the first term in the right-hand side of equation 
vanishes. Moreover from equations (|A1|) we obtain the following first integral "exactly" 



(in the sense that it is a necessary consequence of eq. ||A1|| , with no need for any further 
approximation) : 

h + —fg = constant. (A6) 

On the other hand, if p 3> 1, we can neglect the terms containing p~ 1//2 in the right-hand 
sides of equations ( |A2|) and ( |A3| ) and take either equation in ( |A"l"| ) (but not both) to give 
an "approximate" first integral as (cf. eq. 46 in [Apostolatos et al. 1994| ) 



Xif + X29 = constant. (A7) 



Summing up, we can approximate the system of equations (|A2|) -| 
autonomous differential system: 



df 



15 
64 



X2V, 



dp 

dg 15 

Tp = 64 XlV > 

f p = ~^f-X2 9 )V, 



V 2 



l-f 2 -g 2 -h* + 2fgh, 



with the following 



(A? 



whose integral lines in the (/, g, /i)-space are arcs of parabolas given by the following 
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conditions 

h + \ fg = constant, 
Xif + X2g = constant, 
l-f 2 -g 2 -h 2 + 2fgh>0. 

We can set a new parameter u as 



and write / and g as follows 



f = fo+ , * 2 3 g = go- J 1 «■ (A9) 

Xl + X2 VXi + X2 



Substituting equations ( |A9j ) into equation ( |Aq ) and both equations (|S6|) and (|A9| ) 
repeatedly into equation (|A8|) , we can finally reduce the system of equations (|A8|) to 
quadrature: 



h{u) 



V\u) 



du 



The solution, 



8X1X2 



X1X2 



2(xi + xl) 



u 



xl + xl 
2X1X2 



(X2^o - Xifo) 



(Xl/o + X2go) 2 . l/ , . , x2 , 
1 a" 2 + ^( 2h + fo90) + 

Xi + X2 5 2 

1 



u 



(Xigo - X2/0) 

Vx? + xi 



i 2 



- -(2/i + /ofl-o) 



15 



s Vxl+xi (p -m) = I 



dv 



(A10) 



(All) 



involves an elliptic integral of the first kind, so that u can be expressed in terms of the 
Jacobian elliptic functions ( |Abramowitz fc Stegun 1964j , p. 596). Let us remark that the 
oscillation period is inversely proportional to (xi + xl) 1 ^ 2 an d that the turning points of u, 
i.e., the zeros of V(u), are given by the intersection in the (w, /i)-plane of the parabola 



7 7 (Xi/o ~ X2ffo) 2 X1X2 

8x1X2 + 2(xi + x!) 



u 



X1 + X2 

2X1X2 



(X2^o - Xifo) 
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with the ellipse 



u 



(Xigo ~ X2/0 

Vx? + xl 



+ 5 



1 2 



-(2h + f g ) 
5 



(X1/0 + X2t/o) 2 1, , 



We can also notice that the oscillation amplitude depends on %i and \2 only through their 
relative weight (e.g., the quantity X1/X2 or its reciprocal) and is independent from their 
absolute magnitudes. In fact, the oscillation amplitude is determined only by the right-hand 
side of equation ( |A11|) , in which V(v) is a homogeneous function of xi an d X2 with degree 
zero, as an inspection of equation ( [AIQ ) readily shows. 



Let us finally observe that the approximation condition X\ = X 2 may in practice 
be replaced by the condition \Xi — X 2 \p 1 ^ 2 <C nun{xi,X2} sufficient for the first term in 
equation (|A4j) to be negligible. 



1- Xi 



A. 3.1. Examples. 

'- X2 = X 7^ 0- We have from equation ( |A10| ) : 
h = ho-hil + -Au - u,) 2 , 



V 2 = — 

v 16 



UVA + a 



(A12) 



a + %u 2 ) 



[U — u* 



where u* = (g - fo)/2 1/2 , a = 1 + (3h - f g )/2, A = a 2 + 5V 2 /4. Equation 
( |A11| ) can be integrated now by direct application of one of the formulae on p. 596 
of ( |Abramowitz fc Stegun 1964 ), the selection of which depends on the signs of the 
quantities a — 5ul/8 ± A 1 / 2 . For example, if we suppose that = 0, that is f = g , 
then a + A 1/2 > and a - A 1/2 < 0. The solution is then 



u(p) = -sign(V r 



4(A-a 2 ) 



sd 



A), 



2VA , 



(A13) 



where "sd" is the well-known Jacobian elliptic function, equivalent to the function 
"cn" modulo a translation and a rescaling: sd(v\y) = cn(t> — K(y)\y)/(1 — y) 1 ^ 2 [see 
eq. (16.8.2) in |Abramowitz fc Stegun 1964] ; from now on K(y) is the complete elliptic 
integral of the first kind with parameter y, as defined by eq. (17.3.1) in the last 
reference]. As p varies, u(p) oscillates with a period 



A 



256 #((VA- a)/(2 v / A) 



(A14) 
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The analytic equations ( A12| ), ( |A13|) , and ( |A14j) show the dependence on the initial 
conditions f Q , g , h explicitly. Let us notice that the amplitude of oscillations 
does not depend on x, in this case. A plot of f(p) and V(p), together with the 
numerical solution of equations ( |A2| )-(|A"5]) for comparison, is shown in Figure 6 for 
the choice of po = 75, fo = go = 0.25, h = —0.5, Vo = —0.75, xi — 1> Xi — 1) 
corresponding to Figure 12 of ( [Abramovici et al. 1992|) . Here we have a = 7/32 « 0.22, 
A = 769/1024 w 0.75. In this case we notice that, close to the last stable orbit, 
the agreement between analytical and numerical calculations becomes slightly less 
satisfactory. 

2- X1/X2 — > 0. We have from equations (|A9| ) and ( |A~6| ) that / = fo + u, g = go and 
h = ho — gou/2. Thus, equation ( |A5| ) yields 



V (u) 



b 2 (u - u*f 



where b 2 = 1 + 5^/4, u* = (3g ho - 2f - fo9 2 )/(2b 2 



Vq + b 2 u\. Therefore the 



equation du/dp = — 15x2^(w)/64 is easily integrated and yields 



u(p) =u* (signV r ) T sin 
o 



15 

— bx2(p - Po) + (sign^o)^ 



(A15) 



where simp = bu^/a. As p varies u oscillates with a period given by the following 
expression 

A 



1287T 

15^/1 + 5^/4 X 2 

in perfect agreementfj with formula (B25) of (|Cutler fc 1'lanagan 1994 ). 
3- fo — go — 0. We have from equation ( |A1U| ) that the following relations hold: 

ho + —u , 



(A16) 



(A17) 



V 2 = -T 2 {u 2 + a 2 ){b 2 



u 



where we have defined the following constants: Y = XiXz/iXi + xl) 
a 2 = {2[1 + Qh T + (5 + Ah 2 )! 2 } 1 ' 2 + 2 + Qh T}/{bT 2 ) and b 2 = 



4 There are a few minor typographic errors in Appendix B of (Cutler & Flanagan 1994). Equations 
(B14), (B20) and (B25) should read h 4 = -15(s 2 a 2 , l + L5) / (128 fi), k 2 = 1 - a 2 2 l - (1 - al t )a 2 _ Jvl and 
v 2 = [225/4096+11250^/16384] s 2 /^ 2 respectively. 
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{2[1 + 6h T + (5 + 4/i 2 )r 2 ] 1 / 2 - 2 - 6h T}/(5T 2 ). The solution is 
-sign(Vo)a& / 15 



As p varies u oscillates with a period 



b 



.2 



a 2 + 6 2 



(A18) 



A = - (A19) 

Wxf + xi v /i + 6/ i0 r + (5 + 4/ i 2 )r 2 



The analytic equations ( |A17|) , ( |A18|) , and ( |A19|) show the dependence on xi and Xi 
explicitly. Both the oscillation amplitude and period depend on xi and X2 in this 
case, but the amplitude depends only on F and not on (xi + xV) 1 ^ 2 ^ as previously 
noticed. Let us remark that in the limit Xi ^ X2 we find for equations ( |A18|) and 
( |A19| ) the same expressions given in the example @ by equations ( |A15| ) and ( |A16| ) 



with fo = go = 0, as is easily seen when we notice that, if xi <S X2, then T ~ X1/X2, 



a 



2 4x2/(5Xi)' & 2 ~ 1 — and that lim^o sd(v\y) = sinw. 



We plot in Figure 7 the analytical solutions (solid line) and the numerical solutions 
(dots) for f(p) and V(p), with the chosen values xi = 10~ 3 , X2 = 1, fo = go = ho = 
and po — 250. As we see, the agreement between the analytical and the numerical 
curves is striking. 



A. 4. Very different masses at large separations. 



By differentiating equation ( |A5|) twice with respect to p and making repeated use of 
equations ( |A2| )-(TA~4"D, we obtain d 2 V/dp 2 = F(f,g,h,p) where F is an algebraic function, 



easily calculated. If \Xi — X 2 \p 1 ^ 2 ^> 1, then the leading term in F is much greater than all 
others and the equation for V is approximated by the following one 



d 2 V 
dp 2 




whose solution is 



V(p) = d Ai (-Q 2/3 p) + C 2 Bi (~Q 2/3 p), (A20) 

where Ai and Bi are the Airy functions (fully described in Appendix b of |Landau fc Lifshitz| 
19581 and §10.4 of |Abramowitz fc Stegun~T964T ) , SI = 15\X t - X 2 \/(128r)), C x and C 2 are 
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integration constants such that V(p ) = V and [dV/dp] = V ', namely, 



Co 



v Bi' (-ft 2 / 3 p ) + v ' ft- 2/3 Bi (-ft 2 / 3 



Poj 



Ai (-ft 2 / 3 p ) Bi' (-ft 2 / 3 p ) - Ai' (-ft 2 / 3 p ) Bi (-ft 2 / 3 p )' 

V Ai' (-ft 2 / 3 p ) + V ' ft' 2 / 3 Ai (-ft 2 / 3 p ) 
Ai' (-ft 2 / 3 p ) Bi (-ft 2 / 3 p ) - Ai (-ft 2 / 3 p ) Bi' (-ft 2 / 3 p )' 

) to the lowest order terms in p, the 



We then get, by truncating equations ( |A2| ) 
following approximate solutions: 



M = fo- 
g(p) = 90 + 

h(p) = h + 



15 X2 
128 A x 

128 X 2 

15 X l - X 2 
128 77 



V(y) dy, 



V(y)dy, 



P<> 



(A21) 



y 1/2 V(y)dy, 



where V(p) is given by equation (|A20|) . 

For the evaluation of equations ( |A21| ) we need to know the functions J y Ai (-y)dy 
and J y Bi (—y) dy, which are tabulated in Table 10.12 of ( [Abramowitz fc Stegun 1964| ), or 
asymptotic expansions for them (valid if both p and po are much greater than unity), which 
we adapted from equations (10.4.83) and (10.4.85) of ( |Abramowitz fc Stegun 1964|) and 
give here for reference: 



f Ai (-ft 2/3 y) dy ~ -vr-^ft- 

J PO 

f Bi (-ft 2/3 y) dy ~ +7r- 1/2 ft- 



■7/6 



■7/6 



- 3 / 4 cos f-ft y 3 / 2 + - 
V3 y 4 



y 3/4 sin 



-fty 3 / 2 + - 
3 y 4 



- p 

i 

-I po 

p 

PO 



We also need the following formulae, which we worked out with the help of equations 
(10.4.60) and (10.4.64) in QAbramowitz fc Stegun 19611 ): 

f 

Ipo 



y i/2 Ai (-n 2 / 3 y)dy~- 7 r-V 2 n- 7 / e 



-1/4 

y 1 cos 



-fty 3 / 2 + - 
3 y 4 



/V /2 Bi (-ft 2 / 3 y)dy~ +7 r-V2o- 

. •'PO 



7/6 



-1/4 • 

y 7 sin 



p 

PO 

p 

PO 



If Xi 7^ A2, the condition \X\ — A^lp 1 / 2 ^> 1 for the validity of the approximation 
can always apply at sufficiently large values of p, but the agreement becomes increasingly 
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poorer as p decreases. In fact, at short separations the solution given by equations (|A21 ) 
fails to predict both the period and the amplitude of oscillations to a sufficient degree of 
accuracy. Nevertheless this solution is worth being considered because it is very simple and 
approximately describes the behavior of the system in a large variety of cases. 
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Figure captions. 

Fig. 1. Contour lines cJ/{GM 2 ) = 0.25, 0.50, . . . , 2.50 in the (p, r/)-plane. 

Fig. 2. Evolution of Ltv • o"i from p = 75 to p = 6 for a system with m 2 = 0.13m! with the 
following choice of initial values: (L N • o"i)o — —0.78256, (L^ • <x 2 ) = —1.1177 x 10~ 4 , 
(o-i ■ <t 2 ) = -8.7585 x 10~ 5 , (L N ■ a x x <x 2 ) = -3.6197 x 10" 4 . 

Fig. 3. Behavior of the Kerr parameter as function of p when xi — X2 — 1 {solid lines) and 

(a) mi = m 2 ; f = g = 0.9, h = 0.63 at p = 100. 

(b) m 2 = O.Itoi; /o = —go = 0.9, fo Q = —0.7 at p = 100. 

(c) m 2 = O.lmi; f = -0.42555, g = -0.59746, ho = 0.23651 at p = 180. 

(d) m 2 = O.lSnn ; / = -0.999247282, ^ = -8.44485 x lO" 3 , h = -8.45 x 10~ 3 at 
Po = 75, which is the same choice as for Fig. 2. 

The corresponding spinless cases are shown for comparison (dashed lines). 

Fig. 4. Decrease of the Kerr parameter almost to zero at the end of the inspiral phase is 

obtained only for particular choices of parameters and initial conditions. In this case 
we have m 2 = lOmj, Xi = X2 = 0.351 and fo = —go = —h = 1 at po = 100. 

Fig. 5. (a) Probability density function of cJ/(GM 2 ) at p = 6 computed by the Monte Carlo 
method. Assumptions for the probability distributions of parameters and relative 
orientations of angular momenta are explained in text. A normal distribution with 
same mean and standard deviation is plotted for comparison, (b) Comparison of the 
same curve as in (a) with the probability density function of cJ/(GM 2 ) at p = 6 for 
the spinless case (dashed line), numerically computed. Notice that the figure scale is 
different from (a) to (b). 

Fig. 6. Plots of (a) / = Ljv • §1, (b) V = Ln ■ §1 x s 2 versus separation p, when p = 75; 

fo — 9o — 0.25, ho = —0.5, V = —0.75; Xi — X2 — 1 an d mi = m 2 . Solid lines show 
the analytical solutions, dots show the numerical solutions to equations (p\~2|) - ([A~5D . 

Fig. 7. Plots of (a) / = Ljv • Si, (b) V = Ln • §1 x s 2 versus separation p, when p = 250; 
fo — 9o = h = 0; Vo = —1; Xi — 10~ 3 , Xi — 1 an d m i — m 2- Solid lines show the 
analytical solutions, dots show the numerical solutions to equations (|A2])-([A5|). 
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